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Abstract 

We consider the Stokes eigenvalue problem in a bounded domain of with Dirich- 
let boundary conditions. The aim of this paper is to advance the development of 
high-order terms in the asymptotic expansions of the boundary perturbations of eigen¬ 
values, eigenfunctions and eigenpressures for the Stokes operator caused by small per¬ 
turbations of the boundary. Our derivation is rigorous and proved by layer potential 
techniques. 
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1 Introduction 

The field of eigenvalue problems under shape perturbation has been an active research 
area for several decades. Several related problems belong to Stokes systems, which 
are further subdivided by assumptions on the underlying media and on the Dirichlet 
boundary conditions. The main objective of this paper is to present a schematic way 
to derive high-order asymptotic expansions for both eigenvalues and eigenfunctions for 
the Stokes operator caused by small perturbations of the boundary. The properties of 
eigenvalue problems under shape deformation have been a subject of comprehensive 
studies m [21013] and the area continues to carry great importance [8llini[l2lllll[9|. 
A substantial portion of these investigations discusses the properties of smoothness 
and analyticity of eigenvalues and eigenfunctions with respect to perturbations. 
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Let n C be a bounded open domain with boundary of class We consider 
the following eigenvalue problem for the Stokes system with homogeneous boundary 
conditions: 


—Av + V.p = \v in n 

< V.u = 0 in n (1) 

u = 0 in dQ. 

\ 

Here v = {vi,V 2 ,V 3 ) denotes the velocity field, while the scalar function p is the pres¬ 
sure. 

It is well known that this eigenvalue problem admits a sequence of a no decreasing 
positive eigenvalues 0 < Ai < A 2 < • • • < < • • • tending to infinity as n ^ -|-oo. 

The eigenfunctions {vn}n>i C and the eigenpressures {pn}n>i C L^(H) may 

be taken so that {vn}n>i constitutes an orthonormal basis of 

H{Q) := {v G (L^(fi))^ : V • u = 0 in H, and tt = 0 on 9H}. 

The pressure p is determined up to an additive constant. 

We assume that H has a small and smooth deformation and that the boundary of the 
deformed domain fij is the set of points x defined by 

dfis := {x = X + dp(x)iy(x), x G dfi} (2) 

where iy(x) is the outward normal vector on on dfi and p(x) is a real function in 
that satisfies 

Ilp(2^)llc2(an) < 1- (3) 

Obviously, the domain is of class and the Dirichlet eigenvalue problem for the 
Stokes system can be defined in flj as well. 

In this paper, we derive the asymptotic of eigenvalues, eigenfunctions and the eigen¬ 
pressures solutions to the Stokes system: 

-Avs + V.ps = X5V5 in Qs 

< V.vs = 0 in (4) 

vs = 0 in dQs- 

Here we suppose that the eigenvalue Aq is simple. Then the eigenvalue A^ is simple 
and is near to Aq associated to the normalized eigenfunction vs- 

To the best of our knowledge, this is the first work to rigorously investigate Stokes 
eigenvalue problem in the presence of the perturbation ([ 2 ]) and derive high-order terms 
in the asymptotic expansion of A 5 — Aq and vs — vq when <5^0. However, by the same 
method, one can derive asymptotic formula for the Neumann problem as well. 

Zuazua and Ortega have proved in m the regularity of the eigenvalues and eigen¬ 
functions of the Stokes system with respect to the perturbation parameter, by using 
the Lyapunov-Schmidt method. Their proofs are essentially inspired in the work of J. 
Albert mm for the Laplace operator. Our analysis and uniform asymptotic formulas 
of eigenvalues and eigenfunctions, which are represented by the single-layer potential 
involving the Green function, are considerably different from those in [HE]. 
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This paper is organized as follows. In section 1, we describe the main problem 
in this work. In section 2, we develop a boundary integral formulation for solving the 
eigenvalue problem Q, and we present some preliminary results. In the last section, we 
derive by layer potentials techniques formal asymptotic expansion for both eigenvalues 
and eigenfunctions of the Stokes operator. 


2 Integral equations method 

We now develop a boundary integral formulation for solving the eigenvalue problems 
m and ©• The components of the fundamental Stokes tensor T and those of the 
associated pressure vector P, which determine the fundamental solution (T, P) of the 
Stokes system in are given by (see for instance m) 


r dA rl = l a a 

ij\ ) 471- |3;| ' 47rA^ ^ 

1 Xi 


Pi{x) — 47rbip' 

We recall that the ith row, T,- of T satisfies 


-ATi + V.Pi(x) - ATi = P5{x) in 


[v.r, = 0 

where (e*; i = 1,2, 3) is the orthonormal basis of M^. 


m 


(5) 


( 6 ) 


2.1 The potential theory for the Stokes system 

Let us denote hy ip = {ipi, ip 2 , ifs) a complex vector-valued function in the class C^{dQ). 

The hydrodynamic single layer potential with density ip G is the vector 

function S{X)ip{x) defined by 


S{X)ip{x):= r{X,\x -y\)ip{y) da{y), x € 
Jdn 


(7) 


The pressure term Q corresponding to the single layer potential is the function 
given by 

Qip{x) := [ P{x,y)ip{y) da{y), x€M.^\dn. 

Jdn 

For a cureful stady of these potentials, one can refer to HU, m, ffo]- 

Taking into account the well known properties of Green function T, one obtains the 
result that the pair {Sip, Qip) are smooth functions in each of the domains fl and 
respectively. Also these functions are classical solution to the Stokes system (fTl) . 

The continuity and jump relations of the Stokes surface potentials on the boundary 
dQ are described in the following proposition (see [6] pp. 41-42 or [lU P- 66 ): 

Proposition 2.1 Let ip G C'‘^(9n)^ and let S denotes the surface potential defined in 
&■ Then on the boundary dLl the following continuity and jump relations are satisfied: 

{S{X)ip)\^ = {S{X)ip)\.=S{X)ip 
dT{X, \x - y\) 




du 


'an 


dn{x) 


-ip{y) da{y). 


3 








2.2 The boundary integral formulation of the problem 

In this section, we give a boundary integral formulation in order to solve the eigenvalue 
problems m and (jl]). 

Proposition 2.2 Suppose p satisfies Then, there exists 5o > 0 such that the map 
^ 5 (x) defined by 

^s{x) = X + 6p{x)iy{x) (8) 

is a C'^- diffeomorphism from dQ to dTls for d < Jq- In addition, the following equality 
holds 

det (V'h^s) = 1 + trV(pi^)(5 + ^ [(trV(/9z^))^ — tr(V(pj^)^)](I^ + det {V{pv)) 5^,(9) 

where tr means the trace of a matrix. Moreover, we have trV(p;^) = div{pi'). 

Proof Recall that dQ is and so n{x) is a vector-valued function. Since the 
function p{x) is the map is also C^. A simple calculation yields the equal¬ 

ity Q. Consequently, for <5 small enough the map 'h, 5 (x) is a C^- diffeomorphism from 
dit to dits- □ 

We further denote the reciproque function of Thanks to 

define the operator A 5 (A) as follows: 

A(A)¥^(x) = (5n,(A)v9(^'7^))(T5(x)), ip € {L^{dM)f (10) 

where <So^(A) is the hydrodynamic single layer potential given by ([7]) when we have 

replaced the boundary dTl by dTl^. 

For i,j G {1,2,3}, we can define the j*^-component of the vector-valued function 
A 5 (A) as follows: 

(.4*(A)v)j(x)=((S!!,(A)v).(4'i-‘))(<PjW), V € (L"(8f!))’. (11) 

The j*^-component of the single layer potential (A) is given by 

{Snsi\)p)Ax) := [ Tji{\,\x - y\)pi{y) dasiy), x G j = 1 , 2,3, (12) 

Jaus 

where pi is the T^-component of the vector-valued function p. Using Proposition 12.21 
relations (Illh - (ll2h and the continuity relations given by Proposition l2.ll we obtain for 
x G dTl that 

(A'^(A)(/p) (x) = [ rji{X,\dis{x) - ^ 5 {y)\)det{V^s{y))Ti{y) da{y), j = 1,2,3. 

Jan 

(13) 

Let Aq the operator dehned as in (fTOl) by 

Aofi = 5(A)0, 

where cj) G {L‘^{dTl)Y. Then, we have the following result, which is a slight variation of 
the Lemma 6.1 due to Ammari and Triki ^ on the scalar eigenvalue problem for the 
Laplacian. 
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Proposition 2.3 The operator-valued function ^o('^) : —>■ is 

Fredholm of index zero in C\iM“ . In addition the Dirichlet eigenvalues of the Stokes 
system m are exactly its real zeros. 

From Proposition 12.31 we know that if Aq is an eigenvalue of ([T]) then Aq is a real 
zero of yl,o(A). Moreover, for eo small enough, the function ^(A) is meromorphic in 
where D^q{Xq) the disc of center Aq and radius eo, and Aq is its unique pole 
in Deg. Furthermore we have the Laurent expansion 

Alo'(A) = (A-Ao)-% + Do(A), (14) 

where io : KerAo{Xo) — >• KerAo{Xo), and Ro{X) is a holomorphic function. 

Our main results in this section are summarized in the following theorem. 

Theorem 2.1 Suppose that the eigenvalue Aq of m is with multiplicity 1. Then, there 
exists a positive constant <5o(eo) such that for |(5| < <5o, the operator-valued function 
X !-)■ . 45 (A) has a real zero X{6) in Deg(Ao). This zero is exactly the eigenvalue of 
the perturbed eigenvalue problem and is an analytic function with respect to 6 in 
] — 5o,^o[- R satisfies A(0) = Aq. Moreover, the following assertions hold: 

Aj\X) = {X-X{S))-^m + Rs{X), 

e{S) : Ker{AsiX{6)) ^ Ker{AsiX{6)), 

where Rs{X) is a holomorphic function with respect to (<5, A) €] — 6o,6o[xD^q{Xo). 


3 Asymptotic behaviors 


3.1 High-order terms in the expansion of As 


We now present some basic results related to shape perturbation. To begin, let 
{ti{x),T 2 {x)) be the orthornormal basis of the tangent plan to the surface dH at a 
regular point x. Their cross product is so orthogonal to dll at the point x. By chang¬ 
ing their order, we can assume that ri x r 2 is a vector pointing towards the exterior of 
the surface did. Then dividing it by its length yields the unit normal vector y{x), that 
is 


vq{x) 


Ti{x) X T2{x) 
Tl{x) X T2(x)| ’ 


(16) 


for x € dH. Evidently izq = n, where u was introduced in section 1. 


rf = grad T 5 • n. 

and = 

- grad •' 

Using Proposition 12.21 we find that 



rf = Ti + 5Mti, 

and = 

: T 2 + 6MT2 

where the (3 x 3)- matrix M is given by 



/ di{pvi) 

d2{pvi) 

dsipi^i) \ 

M = di{pv2) 

d2{pv2) 

d3{pn2) 

\ diipvz) 

d2{piX3) 

dsipj^s) / 


(17) 
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with Vi means the i-ih. {i = 1,2,3) component of the vector v. 


For 5 sufficiently small, one can see that the outward unit normal vector to is 
given by 

Tf{x) X T^{x) 


|rf(x) X r2^(x)|’ 

for X € Then, the following asymptotic expansion holds. 


(18) 


Proposition 3.1 Let vq be given by \16i) . Then, the outward unit normal vg{x) to 
dLls 0,1 x, can be expanded uniformly as 


vs{x) = v^^\x)+ X ^ dLl, 


(19) 


n=l 


where the vector-valued functions v^'^^ are uniformly bounded regardless of n. In par¬ 
ticular, for X G dQ 

v^^\x) = v{x), (20) 




1 


In X t2 


[ti X Mt2 + Mri X T2- (vq ■ (ri x Mt2 + Mti x r2))r'o] • ( 21 ) 


Proof. 

Considering the expansions (jl7h for 5 sufficiently small, the relation (I18p becomes 


a + 5b + 5^c 
a + 5b + 5^c 


( 22 ) 


where a, b, and c are vector-valued functions given by 


a = n X T2, 
b = Mti X r2 + n X Mt2, 
c = Mri X Mt2. 

So that, by expanding the quotient (fllh as 5 tends to zero, we get the desired re¬ 
sults. □ 

Next, one can use Proposition 12.21 to get a uniformly convergent expansion for the 
surface element as follows: 

Proposition 3.2 Let y = ^ 5 ( 7 /) where ^ 5 ( 7 /) is given by ^ for y G dLl. Then, the 
following expansion for the surface element das{y) holds uniformly for y G dLl: 

das{y) = det (VT^) da{y) = ( 7 x 0 ( 7 /) ai{y)6 + 772(y)5^ a3{y)5^)da{y), (23) 

where tto = 1, 7 x 1 ( 7 /) = V • (pv), 7 x 2 ( 7 /) = ^ [(trV(pz7))^ - tr(V(pi7)2)], and asiy) = 
det {V{pv)). 
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Set 


X = ^ 5 (x), X G dQ 

y = '^s{y), y e <9^^, 


(24) 

(25) 


and define 

0(a;,y) := ^(^'^(x) - 4'5(?/) - (x - y)) for (5/0. (26) 

Recall that is a vector-valued function on dVi, then 0 is also a function on 

(50)^. Moreover, the following holds. 

Proposition 3.3 The vector-valued function Q{x,y) is on (90)^ and there exists 
a constant C > Q that only depends on hi and p such that 

|0(a;,y)| < C'|x-y|, (27) 

|0 (x,y) • (x-y)l < Clx-yl^ (28) 


for all x,y G d^l. 


Proof. Expression (|26ll shows that 

Q{x,y) = p{x)v{x) - p{y)v{y). (29) 

Since dPl is a surface, there exists a constant C" > 0 such that 

\v{x) -v{y)\ < C'|x-y|, 

\v{x) ■ {x - y)\ < C'\x - y|^, 

for all x,y G dCl. 

The last inequalities and the regularity of p yield the results of the proposition. 


□ 


Now, by using (12611 we get 


X — y = X — y-\-6Q{x,y), {x,y)GdQxdQ, (30) 

and the following results hold. 

Lemma 3.1 Let r > 0 be a fixed real and m > 2 be a fixed integer. The following 
asymptotic expansions 

^i\\x-y\ ^iX\x-y\ 


\x — y\^ \x — yf 


■ (l + 6t[^^ (x, y) + (A; x, y)) , (31) 


n>2 




piX\x-y\_-\ 

_ — .r. 


\x - yp 


Rf )(A;x,y) + 5Rf\x-x,y) + (5"Rf )(A; x, y) (32) 


n>2 


hold uniformly for (A,x,y) G Br{0) x dLl^, where Bj.{0) is a ball in the complex plane 
of center zero and radius r. In addition the functions T^\X;x,y) and Rn^\X]x,y) 
are smooth and bounded uniformly on Br{0) x dLl^. The first coefficients are given by: 

4'' (A;x,y) , 
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and 


X, y) = {x — y) ■ 0(x, y)( — X^Sij — 3iX-^ — + 15iA 


riXj \ e 


iXr 


\ \2r , inOj+rjOi) {ri9j + TjOi) 

+{x-y)-e{x,y){AiX^-l 5 ^-X -5^+7-^-777-^- 7 ^A -6—^+-7--r 

' r [x — y) ■ k3[x,y) r-^ r [x — y) ■ ^[x,y) 

+ 15 [_ 7 !i! 2 + ^ (riOj + rjOi) 


^4 (x — y) • 0 (x, y )7 


+15 


(x-y)- 0 (x,y) 


6 ^J - \-7^ + 


TiXj {riOj + rjOi) 


r2 (x-y)-0(x,y) 


where r = \x — y\, rj = Xj — yj, 9j means the component of Q, and 6 ij means the 
Kronecker index. 


Proof. For m = 2 and x ^ y we have 

\x-y\‘^ = \x-y+ se{x,y)\^ = \x - y\^ (l + 6 cf\x,y) + 6 ‘^c^ 2 \^^y)) ^ 


where 


„( 2 )(' \ _ ^( 2 )/ N _ e{x,y)-e(x,y) 

Cl KX,y) - |x-y |2 ’ ^ 

Proposition 13.11 shows that c^i\x,y) and c!^\x,y) are bounded uniformly on dQ^. 
For m > 2, we have: 

-&{x,y) im 


|x — yP = |x — y + 50(x, y)P = \x — y\'^\l + 5- 


X - y| ' 


where x ^ y. 

Using ([25]) . Proposition 12.21 and Proposition 13.31 we see that the regular vector-valued 
function (x, y) 1 —>■ is well defined on dXl x dXl, and it is independent of 6 . 

Therefore, we can expand 

00 

1^ “ yP = Cq"^\x, y) + ^ (5”c^”*^(x, y) uniformly on dPl x dPt, (33) 

n=l 


where the first coefficients Cg”^^(x,y) = |x — yp, c^^\x,y) = m < x — y,0(x,y) > 
|x - yp-2 and c^™^(x,y) = f |x - yp[|§^P + (m - 2)( < > )^]. 


Moreover, if m is even, then c^\x, y) = 0 for n > m + 1. 

Now combining (|33ll for m = 1 with the well known asymptotic expansion of the 
exponential function, we immediately get 


g*A|a: y\ _ x, y) Uniformly on i?r-(0) x dXl x dQ, (34) 

n=0 
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where the first coefficients Kq{X; x, y) = , and the coefficients Kn{X] x, y) 

can be deduced recursively from Cn^ 

Thanks to relations (f30t) . (f33|) and ([Ml) , we can obtain the desired result given by (f3T]l 
where the coefficients can be deduced easily from and Kn. 


To prove relation (f32]l . we firstly expand 

- 1) + %. (]^)»% - 1) (35) 

- 1) + 8,.S*,(i^)(e‘-'l*-SI - 1). 

To simplify, let us denote r = |x — y|, r* = Xi — yi, r = \x — y\, and fj defined by the 
relation 

(36) 

(37) 


dx/ = 


Then, 


dx,i^) = -3^, and dx^dx.i^) = -3% + 15^. 




Now, by using the following result 

dxj - 1 ) = dxj{iXr)e^^^ = 


(38) 


we get that 


4.8,, (.•5l‘-SI - 1) = .:a[ 8,.(S),«, + = a[(k _ ^ + (39) 


To find the desired result in (l32]) , we may use ((37)) , ([38]) and (1391) to see that the relation 
verifies: 


^i\\x-y\ _ ^ 

dx-dx- —:-r3— 


= riA%-3%-A2^^-7iA^57 + i5^^1e*^-+3%-15^. (40) 

L 7^0 <j^O /y»0 T* * V' ^ 


On the other hand, the components of the vectorial relation ()30p can be given as 
follows 

Xi-yi = Xi-yi + 5ei{x,y). z = 1,2,3 (41) 

where 6 i{x,y) means the component of the vector-valued function Q{x,y). 

Then, by relations (fMl) and (l4T]l we deduce that 

fj=xi-yi = xi-yi + 66i{x,y) = ri + 50i{x,y), z = 1,2,3. (42) 

Using both relations (f30l) and (l4^ . we get the following expansion 

fiXj = -b + 0(2(43) 

where the hrst coefficients: = ViVj, = riOj + TjOi and ' ^j- 
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Now regarding (I33p and using the fact that ^ 0 for each integer m. Then 
can expand 


one 


{x,y) + Y^ {x, y) uniformly , 


(44) 


n=l 


where the first coefficients: (x, y) = (cq™^) ^(x, y), and y) = — (cq™^) ^)(x,y). 

By using (H3]l - (Hi]l . we get that 


-SoS+^MS + E'*”'’*™ 




(45) 


n=2 


where the first coefficients: /fioUU, i/) — as'^^CQ™:) ^(s:, y), and 

d!;i>(n,y) = -4«>(4'"'(4'"')-2)(x.y) +aS«'(eW)-t(n,,j). 

To achieve the proof, we may insert all expansions (p^ . (|4^ (for m G {4,5}), and 
(ll5]) lfor m G {5,6,7} into (P0|) . We get that 


pi\\x—y\ _ 1 °° 

—p:-= Rq^\x, y) + y) + ^ y) 

n=2 


\x - yP 


where the first coefficient: 

a A _ 


<^^(A; X, y) = 


(*i) 




.(5) 


?(*i) 


-1 


= iX5ij{c'Q^) ^-35ij{c^Q^) ^-X^rirjic^) ^ - TiAr^rj (c^,®^) ^ + 15rirj{c^J^) 

+35jj(c[)^^)“^ - 15rirj{cP)~^. 

Using the fact that Cq™^ = |x — y]™, we find that 

/..s _ 1 

4*'^(A;x,y) := Qx.g.y ■ 

Based on (|40]) . we find that 

i?f^'^(A;x,y) = {iXd.jKi^^ - 

+ (iX5ij4^^ - 35ij4^^ - X^44^ - lix44& + 15/3ii?)i^o + . 

So that, by using the fact that Ki = iXQ ■ \ we get 




( 1 ) 


(^-y) 


( 5 )h-l 


o(6)h-l 


i?[''^^(A;x,y) = iA0- ^ ^ ^ (^iX5ij{c4) ^-36ij(4^^) ^-A^rirj(c^'’^) ^-7iArirj(c^°^) 

+15rirj(co'^^)“^^e*^'^o ^ 

-7zA[-rirj(c[)®V^cf^+(u6'j+?’j6'i)(co®V^]+15[-rirj(c[,'^V^cf^+(u6'i+?’i6'i)(c[,^V^])e*^''° ^ 
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+35y(c[,^^) -l5[-rirj{c^J^) + {nej + rjei){c^J^) ^]. 

Using the explicit forms of the coefficients and given above, we immediately 
get the desired expression of Ri^\X; x, y). The other coefficients can be deduced easily 
by the same manner in terms of c^\ □ 

Now using (ITT]l . we obtain for i,j = 1,2,3 that: 

{xi-yi){xj -yj) = go{x,y) + 6gi{x,y) + 6‘^g2{x,y), (46) 


where 

go{x,y) = {xi - yi){xj - yj), gi{x,y) = ei{x,y){xj - yj) + 6j{x,y){xi - ?/*), 

and 

g 2 {x,y) = 9i{x,y) ■ Oj{x,y). 

Now, by using (l23]l . (f33|) and (|46]l . we immediately get 
^X" ^^^ das{y) = (fo{x,y) + 5fi{x,y) + ^5'^fn{x,y)^da{y), (47) 

where the Taylor coefficients r„ can be given explicitly with the aid of 50 ) 51 ) and 52 . 
Next, the following result holds. 

Lemma 3.2 The following uniform expansion holds on R x dO. x dVL: 

00 

{xi - yi){xj - yj)\x - y\e'^^^~ykas{y) = Es{X-,x,y) da{y) = 6^En{X-, x, y) da{y) 

n=0 

(48) 

with 

Eo{X]x,y) = go{x,y)\x - 

and the other coefficients E^ are deduced from those an, Cn \ Kn and gn- 

Proof. From the proof of Lemma l3.ll (for m = 1 ) and from relations (I23j) and (I46h . 
one can get that 

{xi - yi){xj - yj)\x - y\e"^\^~^Uasiy) = (^go{x,y) + 5gi{x,y)+ 

00 00 

^'^g 2 ix,y)^(^'^S'^c^n^^(^'^Kn{X-,x,y)^ ( 0 - 0 ( 5 ) + 50 - 1 ( 5 )+ 

n=0 n=0 

5^ 0 - 2 ( 5 ) + 5^ 0 - 3 ( 5 )^ da{y). 

By collecting terms of equal powers in the above relation, one can deduce easily the 
uniform expansion ()48l) with 

EoiX-,x,y) = 5 o(x, 5 )|x - 5 |e*^l^"^l. 
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□ 


Now, from ([5|) we have 


„iXr piXr 

-Uo^i a.7,a.-1-1 c r . 


Ti,{X, |x -y\) = —[^- +—[i^- 4i^ - SirxiXj] 


AttX ^ 


XiXi 

^4 


.XjXi 


+i;^[ - 3^ + 15i^] + [ - 3^ + 15 :^], 


e 

47rA2 J ' 47rA2 

where r = \x — y\ and Xi = xi — yi. Then, by inserting relations ([3T|) . (1461) . (1471) and 
that given by Lemma 13.21 into above identity, we immediately get the following main 
results 

Proposition 3.4 Let the perturbed boundary dLls defined by Let x and y given 
by (f£^- /fg7|) . and the surface element das{y) given by /flTj) . Then, the components 
hij, 1 < *, j <3 of the fundamental Stokes tensor can be expanded uniformly as: 

Tij{X, \x - y\) das{y) = \x - y\) + (5r|j^(A, \x - y\) (49) 

+ '^6 ^T\f{X,\x-y\)J da{y), y € dn 

n>2 

where the first coefficients: 


45), 


4i 


44 ); 


r^lA, |a: y\) 3^2 ^ ^2 j + 47]-^A2 a^\^o 9o 


47rA 


3* ^ 15 ~ 

^0 + 7 — 


47rA 


47rA^ 


and 


(5). 


r!‘>(A, k - si) := J (4" + fix' -fi- 3ifi) + -1) [rf Ai + rLso] 


r7^(4)ft I rp{4) ~ 1 

■47aT 9.+r. 


El + 


47rA^ 


Ti. 


In Proposition 13.41 the coefficients tq and ri are deduced from ([33|) for m = 5. 

Now, introduce a sequence of components of integral operators {A["‘'^)n>o, defined 
for any ip € L‘^{dLl)^ by: 




'.x)^l rf{x,y):pj{y) daiy), for i,j € {1, 2,3} and n > 0. (50) 


I an 


Using previous results, il is clear that we know explecitely the first terms Af*\ a!'P . 

For any positive integer N, we can by recursive method get the term Then, the 

following theorem holds. 
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Theorem 3.1 Let he the operator defined by Let'ig{x) be the diffeomorphism 
defined as in LemmalKM Let N be a positive integer. There exists a positive constant 
C depending only on N, and \\p\\c^ such that for any (p G L?‘{dLlsY and i,j € {1,2,3}, 
the -component Af defined by satisfies: 

N 

n=l 

where ip = (po^^. 

3.2 Asymptotic expansion of the eigenelements 

To develop asymptotic behaviors of eigenvalues and eigenfunctions with respect to 
the parameter of perturbation 5, we may use the results of Theorem 13.11 Then, the 
following asymptotic expansion related to the operator As appears clearly. 

Proposition 3.5 Suppose that we have all hypothesis of Theorem I3.il Then, the 
operator 3l'^(A) defined by 071) can he expanded uniformly for x G dTl as follows: 

A^{X)ip = A^^\X)ip + 6A^^\\)p: + 5‘^A‘'^\X)ip + ----, as 6^0, (51) 

where ip G Lfi{dLl)^, the -component of the first term is given by 

{Af\x)ip) (x) = [ T\f{x,y)ipj{y)da{y), for i, j e {1,2,3}, 

Jon 

and more generally, the component of the term with order n is given by 

{A^'^'>{X)p>) (x) = [ r^-^\x,y)p>j{y) da{y), fori,j G {1,2,3} and n > 1. 

Jan 

The coefficients are given by [JW - 
Let aj{6) denotes: 

ajiS) = [ (A - Ao)''(Al'')-nA)5AAl''(A)(/pdA, (52) 

JaD,Q 

where ip G L‘^{dFl)^. The functions aj{6) is analytic in a complex neighborhood of 0 
and satisfies: aj{6) = aj{5). The following main result holds. 

Proposition 3.6 Suppose that Aq is an eigenvalue with multiplicity 1 of the eigenvalue 
problem m- Let the operator A^ be defined by ilM . Then, there exists a small positive 
number 5o such that the eigenvalue X{6) is analytic in ] — (5o)<5o[ and satisfy: 

A(3) = Ao + (7Al + ^A„3^ (53) 

n>2 

where the first coefficients are given explicitly by: 

= (A-Ao)[(Al(°))-'(A)Al(i)(A)(Al(°))-nA)5AAl(°)(A)](^dA, 

< A2 = ;^ir [ (A-Ao)[(.4W)-UW(.AW)-i5A.A« + (.4(o))-U(2)(^(o))-iaA^^^^ 

, (>l®)-^(.4(^)(.4(°))-i)^9A>i®]<i^ dX, 
where ip G L‘^{dFl)‘^ 
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Proof. Let Xs be the eigenvalue of the eigenvalue problem ([1]) and let (p € 

Then, if we take the curl of the first equation in ([4]), we see that there exists a function 
w = w{vs) called vorticity associated to vs such that 


Aw + Xsw = 0 in n 
u = 0 in di}. 


(54) 


That is, w is an eigenfunction of the negative Laplacian, but with boundary conditions 
on the velocity vs. 

Next, let us be the stream function for vs given as in Lemma 2.10 of [8]. Then w = Aus 
and Vu ,5 = 0 on d^l. Since us is determined only up to a constant we can then assume 
that = 0 on dQ. 

Thus, Us satisfies the following eigenvalue problem for the Dirichlet biharmonic oper¬ 
ator: 

f AAus + XsAus = 0 in n 

U5 = 0 • 


1 - 


m 


on. 


Note that Temam [16] exploits the similar correspondence between the Stokes problem 
and the biharmonic problem in the proof of the regularity of solutions to the Stokes sys¬ 
tem and to justify several results. Moreover, as pointed out by Ashbaugh in [4], there is 
a similar correspondence between the eigenvalue problems for the Dirichlet Laplacian 
and system (f55|) with the boundary condition V • = 0 replaced by Aus = 0. 

Then, one can exploit this correspondence to use the approach used in |3| to develop 
an asymptotic expansion for the eigenvalue. 

So that, it is well known that there exits a polynomial-valued function 5 i—)• Qs{X) 
of degree 1, analytic in ] — Jq) ^o[ and of the form: 

Qs{X) = X - ai{6) 

such that the perturbation A ,5 — Aq is precisely its zero. For the existence of Qs one can 
follow the general approach used, for example, in |3| for the case of Laplace operator. 
Writing: 

25(^5 — Ao) = 0 . 

Then we have 

— Ao = ai((j). 

Therefore, by (|52h we have 

A 5 - Ao = [ (A - Xo){AYH>^)dxA^{X)ipdX. (56) 

On the other hand, for <5 in a small neighborhood of 0, the following Neumann series 
converges uniformly with respect to A in dDf.^: 


(M'5)-1(A) = (M(°))-^(A) + ^(M(°))-^(A)[(M(°)(A) -M'5(A))(M(°))-^(A)]^ 

So that. 


k=l 


Xs-Xo = ^tr [ (A - Ao)(M(o))-1(A)5aM'5(A)v^ dX+ 
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« oo 

/ (A - Ao) [{A^^Hx) - ^-"(A)) (Al(°))-'(A)] ^dxA^{X)v dX 

k=l 


By using ([HI) , we find that - tr [ (A — Ao)(Al^°^) ^(A)5 aAI‘^(A)(^ dX = 0. This 

2*7r JgD^^ 

result is a direct consequence of the fact that TZo{X) and 5 a^'^(A) are holomorphic in 
the variable A. 


Now we have: 


A5-A0 = ^tr f {X-Xo){A^^Y\X)[{A<^^Hx)-A\X)){A<^^Y\X)]dxA\X)^dX+ 

J dD^Q 

(A-Ao)(Al(°))-nA)[(Al(°)(A)-Al^(A))(^W)-'(A)]'aAAl''(A)^dA+ 

^tr [ (A-Ao) j;(^(°))-'(A)[(^(°)(A) -AI'(A))(AI(°))-1(A)]"5 aAI'(A)v 9 dA. 

JdDef^ ^^3 

Inserting expression (I5ip into above relation, we may get: 

2*7r 

(57) 

+5^dxA^'^^ + ---)<^ dX + ^tr [ (A-Ao)(^®)"^[< 5 >l(^^(^(°))“^+ 


^2^(2)(^(0))-i + ... ] 2 (a^^(o) + + 5‘^dxA^^^ + • • •)(^ dA + • • • 

If we collect the same powers of 5, then we get the desired results. □ 


Define, 

BsiXMx) = (>V(A)(^)(T,-i)(T,(x)), 

where (p € L?‘{dQ) and given by Section (H and W(A) is the operator associated 
to hydrodynamic double layer potential unittiiiiz!. Then, the following main results 
hold. 

Theorem 3.2 Let As be the operator defined by [W\} . and M.s = As + Bs- Let Oq be a 
bounded neighborhood of Ll in M^. Then there exists a constant di > 0 smaller than do 
such that the eigenfunction v{S) corresponding to the eigenvalue, A(d), in (id^(D^))^ n 
H{Lls) can be chosen to depend holomorphically in (x,d) € Oqx] — di,di[. Moreover 
this eigenfunction satisfies the following asymptotic formulae 

u(x;d) = uo(x) + ^u„(x)d”, (58) 

n>l 

where the function vq is the eigenfunction solution of m associated to Xq. The terms 
Vji are computed from the Taylor coefficients of the operator valued function M.s and 
of those of the function a{6) = (ajj(d))i<ij< 3 . 
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Proof. From [iniiii] we deduce that there exists a continuous function which 

is analytic in M^x] — Jq) <5o[ and such that 

n5(x) = 5(A5)(^ + VV’(A5)(/3, xefl (59) 

solves the eigenvalue problem dH). Moreover, the function given by 

satisfies the eigenvalue problem dH) in ^Is with the boundary conditions: U{ 6 )\Qng = 0. 
Here, >^(A(5))(/p(^'-^ <5) = Ms{X)p. 

Now, by dSH), we deduce that v{x; 6 ) = U{ 6 ){x) = Af (A((^))(^('I'“^, 5) is jointly analytic 
with respect to (a:, 6 ) in {||x —d' 5 =o(y)|| < ^ojx] — 6 o, (5o[, where zq is a positive constant. 
The fnnction v{x; 6 ) is jointly analytic in the variables {x,S) G Oqx] — <5o)<5o[- 
We shall now give the asymptotic expansion of the function n(x; 5) when 5 tends to 0. 
Integral equation dZ]) gives us 


v{ 6 ){x) = M{X{ 6 ),\x-^s{y)\)p{y,S)\V'i> 5 {y)\da{y), (60) 

Jan 

where M is the kernel of the operator The perturbed eigenvalue X{5) is in a small 
neighborhood of Aq for small values of 5. Then we have the following Taylor expansion 

M{X{ 6 ), \x - 'k 5 (y)|)|VT 5 (y)| = M(Ao, |x - T(y)|)| VT(y)| +Y,5^ Mk{x,y), 

k>l 

(61) 

which holds uniformly in x € Oq and y € dQ. The analyticity of the function ip{y, 6 ) 
with respect to 5 immediately gives 

(fit, 6 ) = (foiy) + E d^Pkiy), (62) 

k>l 

uniformly in y G dQ. Substituting the last two asymptotics into (|60p we find 

k „ 

u(x; J) = u(x; J = 0) + ^ / ipk-n{y)Mn{x,y)da{y)]. (63) 

fc>l n=l 


□ 


The next result provide us with the asymptotic expansion of the eigenpressures. 

Corollary 3.1 Suppose that we have all hypothesis of Theorem Then the eigen¬ 

pressures ps solution of 0 have the following uniform asymptotic expansion: 

P5{x) =Po{x) + '^pnix)6'^, (64) 

n>l 

where the funetion po is the eigenpressure solution of m associated to Aq- The terms 
Pn are computed from the Taylor coefficients An and Vn = {v\,v‘^,v‘^) as follows: 

Pn{.; Xi; .) = / AuJjdxj + ^Afc / vl^_kdxi, where i = 1,2,3. 

k=o 
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Proof. From system (j3|) we have 


V.ps = Avs + Xsvs 

. Hence, we can expand the function ps in powers of 5 as we have done for A 5 and vs- 
Moreover, we have: 

diPs = Avg + Xgvg, for i = 1,2, 3. (65) 

To get the coefficients of the formula (j64p one can insert both asymptotic expansions 
(f5^ and (f58]l into relation (f65]l . and integrate with the convenable variable. □ 
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